First principle calculations of conductance within plane wave basis set via 
nonorthogonal Wannier-type atomic orbitals 
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We present a plane wave/pseudopotential implementation of the method to calculate electron 
transport properties of nanostructures. The conductance is calculated via the Landauer formula 
within formalism of Green's functions. Nonorthogonal Wannier-type atomic orbitals are obtained 
by the sequential unitary rotations of virtual and occupied Kohn-Sham orbitals, which is followed 
by two-step variational localization. We use these non-orthogonal Wannier type atomic orbitals to 
partition the Kohn-Sham Hamiltonian into electrode-contact-electrode submatrices. The electrode 
parts of the system are modeled by two metal clusters with additional Lorentzian broadening of 
discrete energy levels. We examined our implementation by modeling the transport properties of 
Na atomic wires. Our results indicate that with the appropriate level broadening the small cluster 
model for contacts reproduces odd-even oscillations of the conductance as a function of the nanowire 
length. 



I. INTRODUCTION 

The last decade has witnessed a remarkable miniatur- 
ization of conventional microelectronic devices. If this 
trend is to continue, elements of microelectronic circuits, 
e.g. transistors and contacts, will soon shrink to the 
size of a single molecule. One of the major goals in 
nanotechnology is the construction, measurement and 
modeling of electronic circuits in which molecular sys- 
tems act as conducting elements^. Accurate and repro- 
ducible measurements of current-voltage characteristics 
have been recently reported for atomic wires and sin- 
gle molecules 2 * 3 *^. The experimental progress has been 
accompanied by considerable advances towards density 
functional theory based calculations of transport prop- 
erties of nanostructures. This activity has been largely 
spurred on by development of several electronic structure 
codes for the first principle transport calculations^*^ 

The prerequisite for non-equilibrium Green's function 
calculation of conductance is the partitioning of the 
Kohn-Sham (KS) Hamiltonian into left /right electrodes 
and contact regions. Such partitioning is straightforward 
if one expands the KS wavefunctions as liner combina- 
tions of atomic orbit als*iS*4*^ but it becomes a formidable 
theoretical problem if the plane waves are used for a rep- 
resentation of the KS orbitals. One of the aims of this 
paper is development of a theoretical scheme to the par- 
tition of the KS Hamiltonian within plane wave basis set. 

It has been recently proposed that Wannier func- 
tions can be used to link plane wave electronic structure 
and Green's function transport calculations Wan- 
nier functions are localized in coordinate space and are 
obtained by an unitary transformation of the KS or- 
bitals. The number of Wannier functions is typically 
much smaller than the dimension of Gaussian-type ba- 
sis sets. The additional advantage is that, we treat the 
electrode and the wire at the same level of theory within 
the Wannier functions basis whereas it is not possible if 
one uses Gaussian basis set. The quality of Gaussian- 
type atomic basis set can not be easily controlled, which 



makes the results of calculation sensitive to a particu- 
lar choice of basis functions. In principle, these prob- 
lems do not exist for plane wave basis set, where its 
quality is simply controlled by a single parameter, the 
cutoff energy^! By increasing the cutoff energy, we can 
always get converged results. There are two disadvan- 
tages which makes standard Wannier functions inappli- 
cable for the calculations of transport properties. First, 
Wannier functions are defined for occupied KS orbitals 
and if the same localization scheme is directly used for 
the virtual orbitals it typically leaves them as delocalized 
as they were before the unitary transformation. Sec- 
ond, the centers of these localized Wannier functions 
are not controllable before minimization. To address 
both issues, we have developed the localization tech- 
nique which yields non-orthogonal Wannier-type atomic 
orbitals (NOWAOs) from the plane-wave based KS or- 
bitals. NOWAOs are the maximally localized functions 
defined via the set of unitary transformations of occupied 
and virtual KS orbitals. Our scheme is based upon the 
combination of two localization techniques: Thygesen- 
Hansen-Jacobsen partially occupied Wannier functions 12 
and Mortensen-Parrinello non-orthogonal localization. 13 
Combined together, these two techniques are used to in- 
clude the virtual KS orbitals and to shift the Wannier 
centers from bonds to atoms. 

One additional ingredient, which is necessary for cal- 
culations of conductance, is the Green's function of the 
electrodes. In this paper, we represent the leads by two 
small clusters with additional Lorentzian broadening of 
the energy levels. 

We have implemented the working equations within a 
plane- wave/pseudopotential codei* and we will demon- 
strate the numerical accuracy of implementation for some 
prototypical test examples. The remainder of the paper is 
organized as follows. In Section II, we describe the details 
of our method to calculate the conductance of nanostruc- 
tures, and it application to sodium atomic wires is given 
in Section III. In Section IV, we conclude the paper. 
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and yields the partitioning of the Hamiltonian into 
electrode- wire-electrode sub-matrices: 



H 




(5) 



If the basis set \u> n ) is not orthogonal, the analogous par- 
titioning should be performed for the overlap matrix S 
{Smn — (w TO |w„)). Matrices H and S can be definited if 
indexes n and and m in H nm @ are associated with the 
atomic positions. But this is not the case when periodic 
boundary conditions are employed and the KS orbitals 
are expanded in plane waves: 
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FIG. 1: (a) Schematic representation of an atomic/molecular 
scale junction system for electronic transportation, (b) Sim- 
plified transportation model used in our method, (c) Energy 
levels for junction system containing three-atom sodium wire. 
Density of states obtained by Lorentzian broadening of the 
electrode levels are also shown. See text for details. 



II. METHOD 



where are the expansion coefficients. The plane 
waves exp(iGf) do not have any reference to the atomic 
positions and therefore the partitioning of the KS Hamil- 
tonian into electrode-wire-electrode subspaces can not 
be performed directly within the plane wave represen- 
tation. Several groups have attempted to overcome this 
difficulty by Wannier function representations of the KS 
orbitalsAia 



A. Partitioning of the Hamiltonian 

We begin with the Kohn-Sham equation for the entire 
nanowire junction 

H\il>i) = E^) , (1) 

where H is the Kohn-Sham Hamiltonian and \i[>i) is the 
Kohn-Sham orbital. The next step is the partitioning 
of the system into three parts: two electrodes and wire 
(typically atomic or molecular wire plus parts of the lead) 
as it is shown in Fig. ^ The partitioning of the system 
results into the partitioning of the Hamiltonian to seven 
sub-matrices (left lead, contact, right lead and the lead- 
wire interactions) and it is performed by transforming the 
representation of the KS Hamiltonian from KS orbitals 
to atomically localized basis sets: 

hfc> = X! °*»k> > (2) 

n=l 

where Ui n is a unitary transformation and \w n ) form 
atomically localized complete basis set. The unitary 
transformation Ui n is applied to the KS Hamiltonian 

H = \ipi}Ei(ipi\ = \w n )H nm (w rn \ (3) 

i nrn 



B. Non-orthogonal Wannier- type Atomic Orbitals 

Wannier functions are localized functions which span 
the same space as the eigenstates of a band or a 
group of bands. Traditional Wannier functions are ob- 
tained by transforming Bloch representation to real space 
representation^ in which Bloch vector k is substituted 
by the lattice vector R of the unit cell where the orbital 
is localized. 

\Rn) = f yeWMirfJ 

V Z7r ) JBZ 

= / E^m^-^dk (7) 

{^) Jbz m 

(k) 

where Umn is an arbitrary unitary matrix. The integra- 
tion in eq. J7J) is done in reciprocal space within the whole 
Brillouin zone. The vector k equals zero for disordered 
systems like nanostructures or molecular wire junctions. 
In this case, the Wannier functions is defined via the uni- 
tary transformation of the KS orbitals^ 

\u n ) = ^ Umn\lpm)- (8) 
m 

There are several different schemes to define unitary ma- 
trix U mn and the choice the unitary transformation can 



be tailored to particular applications. Finding the maxi- 
mally localized Wannier functions is pivotal for the par- 
titioning of the Hamiltonian. Although there are sev- 
eral possible ways to define maximally localized Wan- 
nier functions, the method of minimization of the mean 
square spread stands outi 1 ^ 1 ? 

Wannier functions are traditionally constructed only 
from the occupied KS orbitals. Occupied Wannier func- 
tions are located on the chemical bonds, which sometimes 
makes the partition of junction systems difficult. The 
additional complication is that the sum in cq.@ runs 
over all KS orbitals (occupied and virtual) and due to 
the completeness requirement it is necessary to include 
as many virtual KS orbitals in the consideration as pos- 
sible. Therefore, to be used in transport calculations 
Wannier functions should be constructed in such a way 
that (a) they are atomically localized and (b) they in- 
clude both occupied and virtual KS orbitals. The above 
two requirements are interconnected since to get atomic 
centered Wannier functions, we must anyway combine 
unoccupied anti-bonding states with occupied bonding 
states. The extraction of anti-bonding states from the en- 
tire set of virtual orbitals is difficult computational prob- 
lem, since unoccupied anti-bonding states are mixed with 
some scattering states originating from periodic bound- 
ary conditions. 

Following Thygensen, Hansen and Jacobseni^ we begin 
the localization of the KS eigenstates by constructing the 
linear combination of the virtual orbitals 

N-M 

\4>l) = CmlHM+m) , (9) 

m— 1 

where N is the number of KS orbitals and M is the num- 
ber of occupied states. Partially occupied Wannier func- 
tions are written as^ 

M L 

\4>n) = U mn\^m) + ^ U M+l,n\<t>l) , (10) 

m=l 1=1 

with L being the number of unoccupied anti-bonding 
states. The optimal value of L is yet to be determined. 
We minimize the following localization functional 

n = Y,[(k\r 2 \k) - (kWk) 2 ] (11) 

n 

to choose a suitable unitary transformation U mn and to 
obtain the maximally-localized Wannier functions ^Li£ 
Minimization of functional f2 is equivalent to maximiza- 
tion of the functional 

s=EE* 2 - (12) 

n I 

where matrix Z 7 is defined as 

ZL = {4>m\e- idl - ? \4> n ) , (13) 
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FIG. 2: (Color online) Non-orthogonal Wannier-type atomic 
basis functions (NOWAOs) for one gold atom in a simple gold 
junction system. Only half of the junction is shown. 

with Gj and wi being the reciprocal lattice vectors and 
corresponding weights^ For simple orthorhombic su- 
percell, I ranges from 1 to 3, corresponding to x, y, 
and z. In practical implementation, Z 7 is calculated by 
Z 1 = UtZ^U, with Z£ defined as 

(4)™ - ^ m \e- ie '- ? \i>n) ■ (14) 

Analytical gradients of functional S (|12|) arc necessary 
to perform effective maximization. If we write the uni- 
tary matrix at iteration i as Uj = Ui_i exp (—A), then 
the gradients of functional S with respect to A can be 
approximated as 

(d~/dA) tJ =^ Wl [z ji (z; j - zi) - z^Zu - z^)] . 

(15) 

The gradient with respect to the coefficient matrix c; m is 
computed by the following formula 12 

(dE/dc*)^ = (16) 
w/[[Z Udiag(Zt) + Z t Udiag(Z)]Ut] 7V+J . A , +J , 

where diag(7i) is the diagonal part of matrix Z, and U is 
the rotation matrix from KS orbitals to partly occupied 
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Wannier functions with dimension N x (M + L). The or- 
thonormality constraint on matrix c is invoked through a 
set of Lagrange multipliers. The steepest descent method 
is used to maximize S. After the maximization, the un- 
occupied anti-bonding states are obtained from the coef- 
ficient matrix c through Eq.@. 

A final set NOWAOs is computed via the additional 
rotation of partially occupied Wannier functions (f 1 Of) : 



M+L 



ran \ Vm/ 



(17) 



Conductance Formula 



The starting point for the conductance calculations is 
the Landauer formula 19 



2e 2 

G = =^T(E F ), 



(23) 



where T is the transmission function and Ep is the Fermi 
energy of the electrodes. Having obtained the partitioned 
Hamiltonian JJJ, we can compute transmission as the 
trace of Green's function G and coupling matrices / R : 7 



Rotation matrix V nm is defined by minimizing the fol- 
lowing function 13 independently for each n: 



H„ = (u n \p(r- R n )\uj n ) 



(18) 



The weight function p(r) is chosen in such a way that it 
has a minimum at r = to localize each \u) n > around 
R n . Following Mortensen and Parrinelloii we select func- 
tion p(f) as 



,2tt . 
1 - cos(— r Q ) 



Q n can also be written in the matrix form: 

n n = (vtp(")v) rm = 4pWtf„ , 

where v n is the 7ith column of V, and 



R 



(n) 



(4>i\p(f- Rn)\4>j 



(19) 



(20) 



(21) 



with \4> > being the set of partially occupied Wannier 
functions obtained by maximization of S. Matrix element 
p(") are calculated only once and stored for everyn. The 
minimum of Q n is obtained when u n is equal to the nor- 
malized eigenvector corresponding to the smallest eigen- 
value of PW. If we need several NOWAOs for a single 
atomic site, the corresponding number of smallest eigen- 
vectors should be chosen. The number of anti-bonding 
states necessary for the localization (L in ea. (|10f> can now 
be computed by the following formula: 



L = N A N LE - M 



(22) 



where Na is the number of atoms in the system and 
N^E is the number of the lowest eigenstates included 
into Mortensen-Parrinello localization. For example, for 
gold atom, which is typical electrode material in molecu- 
lar electronics, five d-type and one s-type NOWAOs per 
atom are need. In Fig. [5J we plot the generated six 
NOWAOs for a gold atom in a simple gold wire junc- 
tion system. We can clearly see that these six NOWAOs 
reflect the s and d characters of gold atom. 



T(E) = Tr\T L {E)G{E)T R {E)G\E)}. (24) 

The matrices G and T L / R are expressed by the matrix 
blocks of the Hamiltonian H, overlap S, self energy S: 



and 



G(E) = {ES W -H W --£ L --£ R }- 1 (25) 



T L/R (E)=i[V L/R --£{ /R ]. (26) 



The self-energies Sl/^j are defined via the Green's func- 
tion of the left and right electrodes g and the electrode- 
wire interactions: 



V L/R (E) = 
(ES 



WL/R 



H 



WL/RJ&L/ 



g L/R (E)(ESl L/R 



H 



WL/R 



(27) 
)■ 



It is not possible to include the entire leads in the prac- 
tical calculations. The interaction between the wire and 
the infinitely large leads is accounted by the self-energy. 
Different theoretical models have been proposed to ob- 
tain g. The most accurate schemes relay on the sur- 
face Green's function method. Even with some special 
techniques, such as transfer matrhi22i2i and decimation 
technique^, the computation of surface Green's function 
is complex and time-consuming. Since detailed exper- 
imental geometrical structures of leads and a nanowire 
are unclear in most cases, it enables us to use the simpler 
models for the leads. In this paper, the electrode is rep- 
resented by a metal cluster with additional energy level 
broadening^. This model has proven to be very success- 
ful in estimation of the bulk density of states from small 
cluster calculationsi^^ As shown in APPENDIX A, the 
Lorentzian broadening a in the electrode density of state 
is the same as is positive infinitesimal a in the electrode 
Green's function: 

g= [{E + ia^L/R-nL/R]- 1 . (28) 

III. TEST RESULTS 
A. Computational methods 

To illustrate the performance of our method we com- 
puted transport properties of Na nanowire. Nanowires of 
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metal atoms have recently attracted much attention be- 
cause of their fundamental and technological importance. 
In particular, sodium atomic wire has been studied 
both experimentally 26 ! 27 ! 28 and theoreticallyi 29 i 30 i 31 i 32 i 33 
It was found that the conductance of Na wires exhibits 
even-odd oscillation as a function of the number of atoms 
in the wire. The conductance for odd sodium atom num- 
bered wire is close to the quantum of conductance Go 
(2e 2 /h), while the conductance for even sodium atom 
numbered wires is smaller than Go- Different implemen- 
tation and junction models lead to different the values of 
the conductance for even numbered nanowires (0.5 - 0.9 
g o ^ | 29 4 30 i 3i 4 32 J 33 w e use Na atomic nanowire as a proving 
ground for our implementation and we aim to reproduce 
the odd-even oscillation and values of the conductance. 

The calculations were performed using our implemen- 
tation of the formalism presented here in the CPMD 
package 14 . All systems were treated employing peri- 
odic boundary conditions and the KS orbitals were ex- 
panded in plane waves (50 Ry cutoff) at the T point 
of the Brillouin zone. We used local density approxi- 
mation for the exchange and correlation functional and 
Stumpf, Gonze, and Schettler pseudopotentialsM for core 
electrons. The system is simulated by a cluster in a 
large supercell. The size of supercell is chosen in such 
a way that the distance between the nearest atoms in 
the neighboring cells is larger than 8.5A, so that the in- 
teraction between supercell images is negligible. An ex- 
tensive set of the KS virtual orbitals is computed via 
Lanczos diagonalizationS to ensure that all possible un- 
occupied anti-bonding states are included. To speed 
up convergence of the self-consistent iterations, free en- 
ergy functional 3 ^ is used with the electronic temperature 
T = 300 K. Since sodium is a single valent atom, only 
one NOWAO per Na atom is constructed from the KS 
orbitals. 

The whole system is divided into three parts: left elec- 
trode, central wire, and right electrode. The electrode 
part is obtained by cutting a few atoms from Na (001) 
surface. In particular, as shown in the inset of Fig. [3] 
we cut a five-atom cluster, which is composed of square 
four-atom base and apex atom. The geometry of this 
five-atom cluster is fixed to the bulk values. The wire 
part is a single chain of Na atoms, where the distance 
between the atoms is constrained to the nearest neigh- 
bor distance in the bulk system. The distance between 
the electrode part and the wire part is optimized. The 
optimized value d is listed in Table [I] and it shows small 
~ 0.1 — 0.2 A odd-even oscillation as the length of the 
wire varies. We use optimized value of electrode-wire 
separation in all our calculations. 

Let us take a three-atom Na wire as an example to 
illustrate conductance calculations by the method we de- 
scribed in the previous section. 

First, electronic structure of the whole system is calcu- 
lated, and a set of KS orbitals is obtained. The number 
of KS orbitals should be large enough to include all unoc- 
cupied anti-bonding states. Fig|3]shows the conductance 
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FIG. 3: Conductance of three-atom Na atomic wire as a func- 
tion of number of Kohn-Sham orbitals calculated to construct 
the NOWAOs. The inset shows the geometry of this junction 
system. 

as a function of the number of KS orbitals included in 
the localization procedure to define NOWAOs. The de- 
viation of the conductance from the correct value (1 Go) 
vanishes rapidly as the number of virtual KS orbitals is 
increased. Our test calculations illustrate the number of 
the virtual KS orbitals can be adjusted to achieve any 
desired level of accuracy in the conductance calculations. 
The results on FigEl demonstrate that 45 virtual KS or- 
bitals are sufficient to get the converged result for three 
atom Na wire. 

Second, this set of KS orbitals is used to construct 
NOWAOs, with which Hamiltonian matrix H and over- 
lap matrix S are calculated. By solving the correspond- 
ing generalized eigenvalue problem, we get a set of energy 
levels e for the whole junction system, as shown in the 
left column of Fig. The differences between these 

eigenlevels and the Kohn-Sham energies are very small, 
especially for occupied state. At the same time, partition 
of H and S is implemented, and the energy levels of elec- 
trodes, e L and e R , are obtained by solving the generalized 
eigenvalue problem for the corresponding submatrices. 

Third, we introduce Lorentzian broadening of the left 
and right energy levels to continuous density of states, to 
compute the self-energy matrices. The broadening pa- 
rameter a is chosen to be 0.025 Hartree here, and the ef- 
fect of broadening parameter on conductance will be dis- 
cussed later. The transmission probability is computed 
by using section II. B equations. The peak positions of 
the transmission curve can be considered as renormalized 
eigenlevels of the central three-atom Na wire coupled to 
the electrodes (e c in Fig. [TJ. 

In units of Go, the conductance equals to the value of 
the transmission at Fermi energy of the leads. Fermi en- 
ergy for the electrode is not known apriori and it has 
to be computed within the approach. There are two 
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schemes to compute Fermi energy of the electrodes. The 
first one determines the Fermi energy at which the inte- 
grated density of states should be equal to the number 
of electrons in the electrode cluster^ The charge pop- 
ulation of the electrode cluster is calculated by Pl/r — 
^ iei ^(S 1 ^ 2 PS 1 / 2 )ii, where P is the density matrix on 
NOWAOs basis set, which is obtained from the gener- 
alized eigenvalue problem of the whole junction cluster 
system. The other one is more simple, we can just use 
the Fermi energy of the whole junction system. When the 
electrode clusters are very large, both methods become 
equivalent. For the five-atom model used here, we find 
the later method is preferable and it yields 1 Go value of 
the conductance of the three-atom Na wire system. 



B. Oscillation of Conductance 

The most important feature of electronic current flow 
through Na atomic wires is that the conductance oscil- 
lates as a function of the number of atoms in the wire 
We calculate the conductance of Na atomic wires with 
length range from 2 to 5 atoms. The odd-even oscilla- 
tion of the conductance is well reproduced as shown in 
Fig. 0Ji. The conductance for N = 3 and N = 5 is one 
unit conductance Go, while the conductance for N = 2 
and iV = 4 is about 0.7 Go- The transmission curves 
for these wire systems are plotted in Fig. 0Jd. We can 
see that the number of sharp peaks in the transmission 
curves, i.e. the number of resonant states, is equal to 
the number of atoms in the wire. The low temperature 
transport properties are mainly determined by the states 
near the Fermi energy. There is a resonant state at the 
Fermi energy for odd numbered Na wires, but there is 
no such state for even numbered Na wires. This is ex- 
actly the reason of the even-odd oscillation behavior of 
the conductancei^SiS 



C. Eigenchannels 

Since there are multiple peaks in the transmission 
curve, it is interesting to see if sodium wires are single 
channel conductor as we expected for monovalent atomic 
chain. By defining 

t = r|G fl r| (29) 

the transmission can be written as T(E) = Tr(ttt). We 
can decompose the transmission to the combination of 
some eigenchannels diagonalizing the matrix tt^^i For 
sodium wires, we get only one non-zero eigenvalues at 
any energy, which is the manifestation of a single channel 
conduction mechanism. 
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FIG. 4: (a) Conductance oscillation of sodium atomic wire 
with its length, (b) Transmission curves of the sodium atomic 
wires. Fermi energy is set to zero. 

TABLE I: Distance between electrode and wire d (in A), and 
conductance of the Na atomic wire junction systems. The 
length of Na wire N ranges from 2 to 5 atoms. The electrode 
part of the junction system is modeled by one ( El model ) 
or five (E5 model) Na atoms. Conductance for four different 
broadening parameter 0.02, 0.025, 0.03, and 0.04 Hartree are 
listed. See text for more details. 



N 




El Model 






E5 Model 




d 


0.02 


0.025 


0.03 


0.04 


d 


0.02 0.025 


0.03 


0.04 


2 


3.004 


0.33 


0.43 


0.50 


0.57 


3.078 


0.64 0.71 


0.75 


0.77 


3 


3.070 


0.99 


0.99 


0.98 


0.97 


3.270 


0.99 1.00 


1.00 


0.99 


4 


3.024 


0.65 


0.81 


0.92 


1.00 


3.097 0.71 0.72 


0.72 


0.69 


5 


3.042 


0.99 


0.99 


0.99 


0.99 


3.221 


1.00 1.00 


0.98 


0.96 



D. Electrode Cluster Model 

In our method, the electrode is approximated by a 
small cluster with proper level broadening. Therefore it 
is important to understand how the results of our calcu- 
lations are affected by the cluster size and by the choice 
of the broadening. 

The size of the electrode cluster should be big enough, 
so that charge neutrality is maintained for the junc- 
tion system. As it was discussed in the sections III. A 
and III.B the five-atom electrode model is sufficient to 
simulate the transport behavior of sodium wires in our 
method. 38 It is already a very small size comparing to 
typically used models for electrodes, but it is interest- 
ing to see if a smaller cluster still works in our method. 
For this purpose, we also check the simplest one-atom 
model. The results are listed in Table |U Comparing to 
the five-atom model, the oscillation strength of the op- 
timized distance between electrode and wire is smaller. 
We can see that the conductance for one-atom model still 
exhibits the even-odd oscillation, but the conductance of 
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even numbered wires are not very stable anymore. This 
result indicates that the one-atom electrode model pro- 
vides the oversimplified description of the system. 

It is also important to test if the broadening parameter 
strongly effects results of the transport calculations. We 
analyzed the conductances for all wires by varying the 
broadening parameter a from from 0.02 to 0.04 Hatree, 
and list the result in Table [I] We find that for large 
electrode model, the conductance is not very sensitive 
to the variations of a, especially for odd-numbered wire 
system. For one-atom electrode model, the conductance 
changes significantly as a varies. This is because there is 
only one electrode energy level for one-atom model, and 
the shape of the density of states at Ef strongly depends 
on the broadening parameter. 

There is another conductance calculation method for 
single channel nanowirei^S^ It is based on Friedel sum 
rule and also relays on the eigenlevel broadening tech- 
nique to represent continuum of states in the contact. In 
that method, the eigenvalues of the whole junction sys- 
tem e are broadened to obtain density of states, thus the 
broadening parameter should be smaller than the reso- 
nance spacing but larger than single particle level spacing 
for electrode. To satisfy this constraint, the spacing of 
t L and en should be much smaller than that of e c , and 
therefore much larger electrode clusters should be used 
in Friedel sum rule based calculations. In our method, 
however, if the broadening parameter is comparable to 
the resonance spacing, we still get very accurate results. 
Therefore, much smaller metal clusters can be used to 
model electrodes within our implementation. 

IV. CONCLUSIONS 

We have developed and implemented a plane wave 
based method to calculate the conductance of nanos- 
tructures. The fundamental quantity in the present im- 
plementation is NOWAOs which are obtained by the 
multi-step localization of KS orbitals. NOWAOs are used 
to partition KS Hamiltonian to electrode-wire-electrode 
submatrices that is necessary step for the Green's func- 
tion based conductance calculations. The electrode parts 
are modeled by small clusters with proper broadening of 
their eigenlevels and this model is especially suitable for 
systems with very limited knowledge of lead-wire bond- 
ing structure. Transport properties of sodium nanowires 
are studied by this method and the odd-even oscillations 
of the conductance are reproduced. 
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APPENDIX A: LORENTZIAN BROADENING OF 
ELECTRODES LEVELS 

The density of states of the leads is related to imagi- 
nary part of Green's function by 

Img(.E) = -ivD(E) (Al) 

and real part of Green's function can then be obtained 
by Kramers-Koning relation 

Re g (^) = Ip/ +0 °^^. (A2) 

Here P stands for the Cauchy principal value integral. 
For the widely used Gaussian broadening^ this pro- 
cedure is numerically complicated (Romberg integration 
technique) and time-consuming. We use the Lorentzian 
broadening instead and we will show in this appendix 
that the broadening parameter in Lorentzian for the den- 
sity of states is the same as the infinitesimal parameter 
in the non-interacting electrode Green's function. The 
Lorentzian broadening function is written as 

^)-E ^ X+ 2 ■ ( A3 ) 

7T ' (E — tiY + (7 

i v ' 

where q are the eigenlevels for the left or right electrode 
cluster. The Cauchy principal integral can be obtained 
analytically for the Lorentzian density of states and it 
yields the following Green's function matrix for the leads 

or in matrix form 

g(S) = (Jg -H)3 + ^ - H -^ 

= [{E + iajl - H]- 1 . (A5) 
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